Neutrino cross sections and future observations of ultrahigh-energy cosmic rays 
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We show that future detectors of ultrahigh-energy cosmic-ray neutrinos will be able to measure 
neutrino- nucleon cross section, a^N, at energies as high as lO^^GeV or higher. We find that the 
flux of up-going charged leptons per unit surface area produced by neutrino interactions below the 
surface is inversely proportional to (Tvn- This contrasts with the rate of horizontal air showers (HAS) 
due to neutrino interactions in the atmosphere which is proportional to ct^n- Thus, by comparing 
the HAS and up-going air shower (UAS) rates, the neutrino-nucleon cross section can be inferred. 
Taken together, up-going and horizontal rates ensure a healthy total event rate, regardless of the 
value of GvM- 
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Detection of ultrahigh-energy (UHE) neutrinos is one 
of the important challenges of the next generation of cos- 
mic ray detectors. Their discovery will mark the advent 
of UHE neutrino astronomy, allowing the mapping on 
the sky of the most energetic, and most distant, sources 
in the Universe. In addition, detection of UHE neutri- 
nos may help resolve the puzzle of cosmic rays with 
energies beyond the Greisen-Zatsepin-Kuzmin cutoff ||] 
by validating Z-bursts topological defects super- 
heavy relic particles |^ , new strong- interactions |^ , etc. 

Several approved and proposed experiments plan to de- 
tect UHE neutrinos by observation of the nearly horizon- 
tal air showers (HAS) in the Earth's atmosphere result- 
ing from i^-air interactions. The expected rates are pro- 
portional to the neutrino-nucleon cross section. Calcula- 
tions of this cross section at 10^° eV necessarily use an 
extrapolation of parton distribution functions and Stan- 
dard Model (SM) parameters far beyond the reach of ex- 
perimental data. The resulting cross section at 10^° eV 
is ~ 10~'^^cm^ [Q. It has recently been argued that the 
extrapolated neutrino cross section may be too high j^] 
at energies above about 3 x lO^^eV. If the cross section 
is lower, then the event rate for neutrino-induced HAS is 
reduced by the same factor. This reduction would com- 
promise the main detection signal that has been proposed 
for UHE neutrino experiments. On the other hand, the 
extrapolated cross-section may be too low, for it ignores 
possible contributions from new physics. 

We will show, however, that regardless of possible the- 
oretical uncertainties in the cross section, the future ex- 
periments can observe the UHE neutrinos. In fact, a 
smaller cross section would offer a double advantage for 
the planned experiments. First, with a new search stra- 
tegy (described below), the neutrino event rate with a 
small cross section is actually larger than the HAS rate 
with a large cross section. Second, the future detectors 
can also measure the neutrino cross section at energies 
far beyond those achievable in collider experiments. The 
flrst advantage is a boon for neutrino astronomy, while 



the second provides important information for particle 
physics. Here we will take the value of the cross section 
to be a free parameter with a wide range of values. 

This study is motivated in part by a recent analysis of 
upward events by Feng et al. The emphases in the 
two papers are quite different. 

In addition to HAS, proposed cosmic ray experiments 
can also observe up-going air showers (UAS) initiated by 
muon and tau leptons produced by neutrinos interacting 
just below the surface of the Earth, and may possibly ob- 
serve the fluorescence signal from up-going charged muon 
and tau leptons (UCL) themselves. Prior estimates for 
the rate of "earth-skimming" events have used the ex- 
trapolated neutrino cross section A smaller value 
of this cross section reduces the shadowing of UHE neu- 
trinos by the Earth. Therefore, the neutrino angles with 
respect to horizon need not be so "skimming." More im- 
portantly, the expected rate of UCL and UAS may (i) be 
larger, and (ii) depend on the cross section. 

Indeed, a lower cross section increases the UCL rate 
per surface area as cr~^ as long as the neutrino absorption 
mean free path (MFP) in Earth is small in comparison 
with the Earth's radius, i?®; i.e., for a,yN ?J 2 x 10~^'^cm^ 
the UCL event rate is proportional to F^joi/jsi, as shown 
in Fig. 1. This inverse dependence of the UCL rate on 
the cross section is to be contrasted with the rate of HAS 
events resulting from v-wn interactions, which decreases 
as (Ti/AT decreases. Furthermore, we find that for small 
but possible values of the cross section, the UAS rate 
can exceed the HAS rate by several orders of magnitude, 
as displayed in Fig. 2. 

There is a firm prediction for a flux Fi, of GZK neu- 
trinos in the energy range 10^^ to lO^*' eV, based on 
the observed flux of UHECR protons at the GZK limit. 
This flux is expected to peak in the decade 10^^ to 
10^* eV for uniformly-distributed proton sources, and 
around 10^^ eV for "local" sources within ^ 50 Mpc of 
earth |0 . If the GZK flux is the dominant source of UHE 
neutrinos, then one can use the predicted GZK flux value 
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FIG. 1. The ratio Tt of the upward going r flux 
to the incident tau neutrino flux F^^ as a function of 
= ^i^/R(B = l/(o"i'jv^^®)i from eq. (5), with fixed 
\t/R(b = 3.5 X 10~^, appropriate for events initiated by 
~ 10^" eV neutrinos. Here n is the mean nucleon number 
density. Assuming a monotonic cross section dependence on 
^/s, the value of ^ is limited from above by the HERA mea- 
surements, as shown by a vertical dashed line. 

and either the inverse relation between the UCL rate and 
(shown in Fig. 1), or the hnear relation between the 
HAS rate and a^N to infer from future data the value 
of the neutrino cross section at ultrahigh energies. In 
addition, a comparison of the UAS and HAS rates, or 
a measurement of the angular distribution of UCL/UAS 
events may allow an experiment to determine a^N inde- 
pendent of Fi,, as discussed below. More promising for 
determination are possible neutrino fluxes at and 
above lO^o eV given, e.g., in ||,[2). 

Let us now estimate the rate of upward showers. UHE 
neutrinos are expected to arise from pion and subsequent 
muon decay. The initial flavor composition is there- 
fore and Ve with a ratio 2:1. These flavors oscillate 
and eventually decohere during their Hubble-time jour- 
ney. The resulting neutrino state includes a Vr fraction 
1 \U^j? |f/r,f + 1 We,? \Urj?, whcrc are the 
mixing elements relating the neutrino mass and the fla- 
vor bases. If \Ur3,\ — jC^^sl are large, as inferred from 
the Super-Kamiokande data, then the oscillations nearly 
equalize the number of UHE neutrinos of each flavor . 
In particular, F^^ w ^F^, is expected. The energy-loss 
MFPs, A,- and A^, for taus and muons to lose a decade 
in energy are 11 km and 1.5 km, respectively, in sur- 
face rock with density psr = 2.65 g/cm'^, and 2.65 times 
longer for lepton trajectories passing through ocean wa- 
ter. Tau and muon decay MFPs are long above 10^* eV: 
cTr = 490 {Er/lO^'^eV) km, and cr^ ~ 10^ cTr for the 
same lepton energy. Because the energy-loss MFP for a 
T produced in rock or water is much longer than that of a 
muon, the produced taus have a much higher probability 
to emerge from the Earth and to produce an atmospheric 
shower. Thus, the dominant primary for initiation of 
UAS events is the tau neutrino. In what follows we focus 
on tau neutrinos incident at 10^° eV. 

Let us consider an incident tau neutrino whose trajec- 
tory cuts a chord of length I in the Earth. The probability 



for this neutrino to reach a distance x is Pu{x) = e"^/^", 
where A"""^ = g^n P (the conversion from matter density 
to number density via TV^/gm is implicit). The proba- 
bility to produce a tau lepton in the interval dx is 
The produced r carries typically 80% of the parent neu- 
trino energy; we approximate this as 100%. Then the 
T produced at point x emerges from the surface with 
energy Er = E^, e~('~^)/-^^. The probability of a t pro- 
duced at point X to emerge with sufficient energy i?th to 
produce an observable shower can be approximated as 
^'r^uAS = B(A^ +X-1), with A^ = \yi{E^ / Ea,); 

(3t « 0.8 X lO^^cm^/g ||l^ is the exponential energy- 
attenuation coefficient. For taus propagating through 
rock, one can take Ar ~ 22 km for E^ ~ lO^^eV and 
i?th ^ lO^^eV, while for taus propagating through ocean 
At is 2.65 times larger. The rates in Fig. 2 are shown 
for UAS over land. An UAS over an ocean, where part 
or even all of the tau path is in water |l^], is best ad- 
dressed with a simulation. The threshold Et\i depends 
on details of detector sensitivity and aperture, as well 
as the assumed neutrino spectrum; these can suitably be 
incorporated in a threshold parameter. 

Taking the product of these conditional probabilities 
and integrating over the interaction site x we get the 
probability for a tau neutrino incident along a chord of 
length I to produce an UCL: 

P.^^ril) = t ^ e"-/^^ = (e^^A-- - 1) e~'/^- . (1) 

The emerging tau decays in the atmosphere with prob- 
ability Pd — 1 — exp{—2R^H/cTTl), where H ^ 10 km 
parametrizes the height of the atmosphere. Thus, the 
probability for a tau-neutrino to produce an UAS is 

P..^UAS(0 = (1 - e-2«eff/-.;) p^^^^i^i) . (2) 

Next we calculate the probability for an incident neu- 
trino trajectory to have chord length I. Due to isotropy 
of the neutrino flux, it is enough to consider incident neu- 
trinos with parallel trajectories. The length I and impact 
parameter h with respect to the Earth's center of a chord 
are related by ^^/4 -t- ft,^ — i?^. The fraction of neutrinos 
with chord lengths in the interval {I, I + dl} is therefore 

2n h dh I 

Pc^or^m = ^ ^dl. (3) 

To get an event rate probability from the incident neu- 
trino flux, there are two further geometric factors to be 
included: the solid angle factor tt for a planar detector 
with hemispherical sky-coverage, and the tangential sur- 
face area A of the detector [|l6| . 

Putting all probabilities together, we arrive at the rate 
of UCL and UAS events: 

R.^vKs.^F-r^A __P^^_^^(U^g)(Z). (4) 
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FIG. 2. The air shower probabihty per incident tan neu- 
trino RuAS I Fv^ ttA as a function of the neutrino cross section 
(eq.(6)). The incident neutrino energy is 10^" eV and the as- 
sumed energy threshold for detection of UAS is Eth = 10^* eV 
for curve 1 and 10^^ eV for curve 2. 



The double integral in eq. (Q) is easily done analytically 
for the UCL case, and for the UAS case when the angle 
above the horizon satisfies 9 ^ (i?T-/10^^eV) deg, such 
that 2RqH/cTt-1 <C 1. These results in the limit X^, ^ 
At- are 



Rr = F.^TtA -^ri. 
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where = Xt/R® and ^i, — \u/R®, and 



R^,s=F.^^A^ (1 
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For neutrino trajectories through the Earth's mantle, 
= 0.66/cr33, where (T33 is the neutrino cross section 
in imits of 10~^'^cm^. 

The tau flux scales as oc cf^n for X^, ^ i?^ 

(i.e. small ct^at), because the large neutrino MFP ex- 
ceeds the Earth's diameter, making the interactions rarer 
for increasing Xu- For Ar <C Ai/ <C i?e (i.e. cfvn ~ 
2 X 10^'^'^cm^), the probability scales as X^, oc (J^^- The 
rise with increasing MFP is attributable to shrinkage of 
the Earth's "shadow" and the consequent increase in the 
target volume. The ratio r,- = Rt/F^ttA is shown in 
Fig. 1. We note that the neutrino cross section may 
be determinable from the angular distribution of UCL 
events alone (if they can be observed), independent of 
the neutrino flux. One expects the angular distribution 
of UCL to peak near cos^poak ~ Xi,/2R^, which implies 

a^N ^ {'^{P) R® COS^pcak)"^ ■ 

The UAS rate scales differently from the UCL rate, due 
to the dependence on path length in the atmosphere. In 
Fig. |2| we show that the the number of expected UAS 
events per incoming neutrino as a function of the neu- 
trino cross section (using the exact integral expression). 
The cross section is bounded from below by the value 
~ 2 X 10"34cm2 measured at HERA at = 314 GeV, 
corresponding to a laboratory energy = 5.2 x 10^^ eV. 
For comparison, we also show the number of expected 
HAS events per neutrino that crosses a 250 km field of 
view, up to an altitude of 15 km. It is clear that for the 



smaller values of the cross section, UAS events will out- 
number HAS events, while the larger values of the cross 
section favor the HAS events. The ratio of HAS to UAS 
rates may provide a good measure of the cross section. 

We give some examples of the UAS event rates ex- 
pected from a smaller neutrino cross section at 10^" eV. 
Let us choose cr^jv = 10~^'^cm^, for example. Tak- 
ing the mantle density of pm = 4.0(7/cm'^ and i?® — 
6.37 X lO^cm, one gets = 0.65. Reference to Fig. 1 
then shows that the — s- t conversion probability is 
Tt — 0.1% for land events with Er > 10^^ eV, initi- 
ated by a ~ 10^" eV primary neutrino. Including the 
probability for a tau to decay in the atmosphere, the 
Vr UAS probability is 4 x 10~^ (7 x 10"^) for a shower- 
energy threshold i^th = 10^* eV (10^^ eV), according to 
Fig. 2. EUSO and OWL have shower-energy thresholds 
~ lO^^eV corresponding to curve 2 in Fig. 2. They have 
apertures ~ 6 x lO'^km^ and 3 x lO^km^, respectively, for 
a wide angular-range of UAS. These detectors should ob- 
serve F20 and 7F20 UAS events per year, respectively (not 
including duty cycle); here ^20 is the incident neutrino 
flux at and above 10^° eV in units of km~^sr~^yr~^, one- 
third of which are j^'t-'s. Including showers from taus orig- 
inating outside the field of view, and direct tau events, 
increases these rates. The rates may be further increased 
in space-based detectors by tilting towards the horizon 
so as to maximize the acceptance for events with smaller 
chord lengths (where the neutrino attenuation is less and 
the field of view is greater, but the energy threshold is 
higher) and to allow more atmospheric path length for 
tau decay. The rates will also increase if i?th can be re- 
duced; reducing E'th to 10^^ eV or less would allow the 
GZK fiux to be observed and measured. Although the 
spectrum of UHE neutrinos is not known, a variety of as- 
sumed spectra can be encompassed by our parametriza- 
tion, and our conclusions are not sensitive to spectral 
details. 

Finally, we comment on two important inferences. 
First, the reported bounds on the UHE neutrino flux due 
to the non-observation of neutrino-initiated HAS ||l^ and 
of radio signals produced by neutrino interactions near 
the surface of the moon ||l8| are weaker if the cross sec- 
tion is smaller. Concerning the lunar radio bound, the 
lunar radius is 1740 km, about 3.5 times smaller than 
that of Earth, and the density of the Moon is about the 
same as the Earth's surface density. Thus, — A^/i? 
for the moon is 3.6/(733, which is about 5.5 times the 
value for earthly neutrinos. Consequently, the range of 
S^i, for the Moon allowed by our ignorance of the true 
neutrino cross section is very large, and the true neu- 
trino flux limit from lunar radio could be different from 
that previously reported. The second inference has to 
do with the predicted angle- independence for upgoing r- 
neutrinos at '-^ 10^^ eV [Q. For a smaller cross section, 
a harder spectrum of unattenuated j^t-'s above 10^* eV 
and a larger angle dependence may be expected. 



3 



To conclude, the overall prospects for UHE neutrino 
astronomy are not diminished by the theoretical uncer- 
tainties in value of the neutrino-nucleon cross section. 
UHECR neutrinos can be observed at a healthy rate for 
any allowed value of the cross section. Furthermore, 
future neutrino cosmic-ray experiments can determine 
the neutrino-nucleon cross section at energies as high as 
10^^ GeV, or higher, by comparing the rates of UAS with 
those of HAS; or by measuring the angular distribution 
of UCL events. 
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